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ABSTRACT 

We revisit the dynamics of Prometheus and Pandora, two small moons flanking 
Saturn's F ring. Departures of their orbits from freely precessing ellipses result from 
mutual interactions via their 121:118 mean motion resonance. Motions are chaotic 
I because the resonance is split into four overlapping components. Orbital longitudes 

I were observed to drift away from predictions based on Voyager ephemerides. A sudden 

' jump in mean motions took place close to the time at which the orbits' apses were 

antialigned in 2000. Numerical integrations reproduce both the longitude drifts and the 
. jumps. The latter have been attributed to the greater strength of interactions near apse 

antialignment (every 6.2 years), and it has been assumed that this drift-jump behavior 



in 
o 
o 

(N 

> 
O 



X 



will continue indefinitely. 

We re-examine the dynamics of the Prometheus-Pandora system by analogy with 



O 

I that of a nearly adiabatic, parametric pendulum. In terms of this analogy, the current 

^ ' value of the action of the satellite system is close to its maximum in the chaotic zone. 

Consequently, at present, the two separatrix crossings per precessional cycle occur close 
to apse antialignment. In this state libration only occurs when the potential's amplitude 
is nearly maximal, and the "jumps" in mean motion arise during the short intervals of 
I libration that separate long stretches of circulation. Because chaotic systems explore 

the entire region of phase space available to them, we expect that at other times the 
Prometheus-Pandora system would be found in states of medium or low action. In a 
low action state it would spend most of the time in libration, and separatrix crossings 
would occur near apse alignment. We predict that transitions between these different 
states can happen in as little as a decade. Therefore, it is incorrect to assume that 
sudden changes in the orbits only happen near apse antialignment. 
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1. Introduction 

Discovered by the Voyager spacecraft in 1980 and 1981, Prometheus and Pandora are two small 
moons of Saturn. Properties of their orbits are summarized in Table 1. A series of observations 
starting in 1995 found that the orbital longitude of each satellite deviated by about 20 degrees 
from its value as predicted by the Voyager ephemeris (Bosh & Rivkin, 1996; Nicholson et al., 1996; 
McGhee et al., 2001). Further, each satellite's mean motion underwent an abrupt change in 2000 
(French et al., 2002), seen as "kinks" in their mean longitudes. Goldreich & Rappaport (2003a; 
GR03a) pinpointed chaos due to the satellites' mutual gravitational interactions as the cause of 
these discrepancies. The system was found to have a Lyapunov exponent ~0.3 yr""^. Jumps in 
the mean motion were attributed to the stronger interactions that occur when the orbits' apses are 
antialigned. Orbital integrations in GR03a reproduce the observed gradual drifts away from the 
Voyager-hased predictions. Kinks in the mean longitudes are apparent every 6.2 years, at times of 
apse antialignment. Goldreich &: Rappaport (2003b; GR03b) captures the essential dynamics of the 
system by including only interactions due to the 121:118 mean motion resonance. More complete 
numerical simulations of the system, which include the influences of other Saturnian satellites, have 
been performed by Cooper &; Murray (2004); Jacobson &; French (2004); and Renner, Sicardy & 
French (2005). These confirm that the simplified dynamics in GR03b is sufficient to describe the 
chaotic motions of Prometheus and Pandora. 

Physical explanations for the mean motion jumps have not advanced beyond those proposed 
by Goldreich &: Rappaport: all of the above papers cite enhanced interactions at apse antialignment 
as the reason for the kinks in the mean longitudes observed in 2000. They do nothing to dispel 
the expectation that the mean longitudes will continue to display the same drift-kink behavior 
indefinitely into the future. A practical consequence of this belief is that orbits fitted to Cassini 
data avoid times around apse antialignment and assume freely precessing ellipses between these 
times (Porco et al. 2005). 

We show that a proper understanding of the dynamics of the Prometheus-Pandora system is 
more subtle than previously recognized. In doing so we exploit an analogy with the dynamics of a 
parametric pendulum to reinterpret numerical integrations presented in GR03b. We confirm that 
scparatrix crossings are the cause of chaos and that their rate determines the magnitude of the 
Lyapunov exponent. Our focus is on the precessional phase at which separatrix crossings occur. 
Currently the two that take place during each 6.2 yr precessional cycle occur near apse antialign- 
ment and the kinks arise during intervals of libration separated by long stretches of circulation. 
More generally, separatrix crossings can occTir at any precessional phase, with large changes in 
these phases predicted on timescales of one to two decades. When crossings take place near apse 
alignment, the system remains chaotic but displays qualitatively different features than at present. 

This paper is arranged as follows. In §2 we demonstrate the similarity between the dynamics of 
the Prometheus-Pandora interaction and that of a parametric pendulum. We study the behavior of 
an adiabatic parametric pendulum in §3. In §4 we apply these findings to the Prometheus-Pandora 
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system, and discuss both the imphcations and the limitations of our analogy. We conclude in §5. 



2. Similarity of equations 

GR03b showed that the chaotic motions of Prometheus and Pandora arise from their 121:118 
mean motion resonance. Differential precession of the two eccentric orbits splits the resonance into 
four discrete components whose arguments are tpq = tp — 6q, with 

^ = 121A - 118A', (1) 

where A and A' are respectively the mean longitudes of Prometheus and Pandora, and 

Si = 3w, S2 = 2w + w' , 62 = tu + 2w' , 64 = Szu'. (2) 

GROSb note that because of the rapid precession caused by Saturn's oblateness, interactions between 
the satellites produce negligible effects on their apsidal angles and orbital eccentricities. Therefore 
it is adequate to only treat changes in mean motions, or equivalently in the angle ip. 

The equation of motion for ij; reads (GROSb) 



cP'ip /soTn ^ , am' 
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J2CqSmi^-Sg), (3) 

-> q=l 



where M = 6 x 10^^ g is the mass of Saturn, and where the values of terms in Eq. 3 are given in 
Tables 1 and 2. 



Equation 3 can be written as 



-uj^^A{t)sm[i;-^{t)], (4) 



at time t, where loq — 3.8 yr ^ and \A\ < 1. Defining = tp — (j) and thus transforming to the 
frame drifting with the potential at 0, we obtain 



= -Wo^(t) sin * + (?!). (5) 



Provided (j) <€i cOqA, we can drop the last term, leaving 



= -u;5^(i)sin*, (6) 



precisely the equation of motion of a parametric pendulum. The potential in which moves is 
given by 

V{'^,t) = -u^Ait) cos (7) 



-4- 



Plots of ^, <p and are displayed in Figure 1.^ 

The timescale for variation of \V{t)\ is 2Tr/\'UJ — w'\ 6.2 yrs, and the fractional variation in 

- 1/2 

its amplitude is of order unity. The typical oscillation period of is w ~ 27r/|y| ~ 2 yrs, so the 
slowness parameter e ~ 1 x 2/6.2 ~ 0.3. Thus the Prometheus-Pandora system is only marginally 
adiabatic. Rather than examine it directly, it proves instructive to consider a parametric pendulum 
for which adiabaticity is robust. For this we choose 

A{t) = l + acos{Kt), (8) 

where a < 1 and k -C wq, so that for small oscillations the system is adiabatic with slowness 
parameter 

e ^ an/ujo <^ 1. (9) 
Without loss of generality we can set uq = 1. This potential is illustrated in Figure 2. 

2.0.1. Integration 

All integrations in this paper are done in Mathematica 5 using a Runge-Kutta routine with 
fixed step size. Care is taken to maintain numerical accuracy at a level such that differences in 
initial conditions determine the rates at which neighboring trajectories diverge. 

3. Behavior of the simplified parametric pendulum 

A simple pendulum can exhibit motion of two kinds:^ libration for E < \V\ in which is 
trapped in a valley, and circulation for E > \V\ in which ip passes over the hills. Between these 
two regimes lies the separatrix, a singular trajectory of infinite period with energy 

E = ip'^/2- Acosi/j = A. (10) 

As A{t) varies, the motion of a parametric pendulum may switch between these two regimes. 
Transitions from large amplitude librations to circulations may occur as the potential shallows. 
Likewise, transitions in the opposite direction can take place as the potential steepens. No matter 
how slowly the potential varies, adiabatic invariance is violated during separatrix crossings. 



^Although 0, we approximate 4> to be constant at —50.7 yr~^ when we calculate the action for the Prometheus- 
Pandora system. The pendulum analogy takes us a long way toward understanding the behavior of the system. We 

assess its validity in §4.2. 

^for a simple pendulum, A is independent of time. 
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For an adiabatic system the action, J{t), remains nearly constant as A{t) varies. We calculate 
J{t) by freezing A(t) and integrating over a period;^ 

J{t) = jpdq = j^dil: = j \/2[E{t) - V{il)')\ d^;'. (11) 

In a similar manner, we obtain the period from P{t) = ^ dq/q. On the separatrix, 

Jsep(i) = mtf'\ (12) 

Figure 3 illustrates some aspects of the behavior of a parametric pendulum characterized by k = 
0.01, a = 0.2, i.e., slowness parameter e ~ 0.002. We see that J{t) maintains a nearly constant 
value as both the action on the separatrix and the period undergo large variations. Jumps in action 
at separatrix crossings are displayed at higher resolution in Figure 4 along with the evolution of 



3.1. Separation in phase space and Lyapunov exponent 

The rate at which neighboring trajectories in phase space separate is of practical interest 
because it limits our ability to make predictions about the future. For chaotic systems, small 
errors in initial conditions amplify exponentially, invalidating long-term predictions. To track the 
separation of two neighboring trajectories, we Taylor expand the equation of motion (Eq. 6) to 
first order in A'i/^: 



The phase space separation, 



S{t) 



=- A{t)u^ All; COS (13) 
2 /dAV'V 



AV'" 



\ dt J 



(14) 



is calculated by the simultaneous integration of Eqs. 6 and 13. The Lyapunov exponent, £, is the 
limit as t ^ oo of {In S)/t. From Figure 5 we find £ ~ 0.0022, which is close to the rate of separatrix 
crossings. 

The Lyapunov exponent measures the average exponentiation rate of the divergence of neigh- 
boring trajectories. Figure 5a shows how this divergence proceeds through a few separatrix cross- 
ings. At each crossing the separation S between trajectories undergoes a sudden jump, correspond- 
ing to the differential phase delay introduced there between neighboring trajectories. Between 
crossings S{t) displays an oscillation superposed on a linear trend. The oscillation reflects the 
periodic motion of the pendulum, and the linear variation is due to the period difference between 



^When the system is in libration, we integrate over only half a period of the motion, because otherwise the action 
changes by a factor of 2 across the separatrix at fixed energy. 
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neighboring trajectories. Both the jumps in S at separatrix crossings and the hnear variations 
between them produce an increase of S on average. Since on average S and its slope change by a 
constant fractional amount at separatrix crossings, S{t) grows exponentially. 



3.2. Excursions in the chaotic zone 

Only trajectories that come close to the separatrix are chaotic. The range in action over which 
there is chaotic behavior is 

AJ,^c,8{Alll-Alll)^8a. (15) 

A system that starts with action outside this zone never crosses the separatrix and undergoes 
regular motion. A system that starts within it eventually explores the entire range of chaotic 
actions. Coverage of phase space is uniform in the chaotic region. 



Adiabaticity is violated along trajectories that cross a separatrix and the action jumps by 

(16) 



AJ= -44ln 
A 



2 cos — sm — 
2 2 



where tps is the phase of the pendulum at which E = Egcp (Timofccv 1978). The validity of Eq. 
(16) requires that A be constant while the trajectory crosses the separatrix. However, A varies on 
the timescale St = \A/A\ = \ tan Kt/K\ which vanishes at the extrema of A. Thus Eq. (16) does not 
apply within boundary layers of width SJ at the top and bottom of the chaotic zone. We estimate 
SJ CaK^, where C is a (large) dimensionless number. 

Figure 6 displays the quasi-random wandering ( "adiabatic chaos" , Neishtadt 1991) of the action 
of a model system with a = 0.2, k, = 0.1 during approximately 5000 separatrix crossings. "Sticking" 
near the top and bottom of the chaotic zone is a consequence of the smaller size of the action jumps 
in these regions as predicted by Eq. (16); A/ A = at the extrema of A{t). 

Consider the distribution of jump sizes. When crossing a separatrix, the system trajectory 
passes over a potential peak. Smaller margins of clearance correspond to slower passages and 
result in larger jumps in the action. The term sin('i/^s/2) is proportional to the difference between 
the system's energy and the separatrix energy when the system makes its last passage through a 
potential trough prior to crossing a separatrix. We construct the probability distribution P(A J) of 
jump sizes under the assumption that sin('05/2) is uniformly distributed between ±1. Over many 
separatrix crossings the action performs a quasi-random walk in which the diffusion constant is 
proportional to f d{AJ)P{AJ)(AJ)'^ . From the plot of the integrand in Figure 7, we see that jumps 
larger than 8A/A, which comprise about 4% by number, contribute about half of the diffusivity. 

Integrating over jump size, we obtain a diffusion coefficient: 

0.82 — -, (17) 
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where tt/k is the mean time between separatrix crossings. Taking values appropriate to the center 
of the chaotic zone, we estimate the zone crossing time 

tc~^~ 15,000, (18) 

in reasonable agreement with that observed in Figure 6. 

The pendulum mostly librates or mostly circulates, according to whether its action is close to 
the minimum or maximum value in the chaotic zone. This is illustrated in Figure 8. As expected, 
the system samples each type of behavior. 



4. Prometheus and Pandora 

Now we return to the Prometheus-Pandora system and detail similarities between its behavior 
and that of a parametric pendulum. Although our interpretation of the satellites' motions rests 
on the pendulum analogy, we integrate the full equation of motion (Eq. 5) for ip(t) and the linear 
equation for the tangent vector Aip{t) derived from it. Because the system is chaotic, it explores 
the entire accessible region of phase space for arbitrary initial data. We find it most convenient to 
adopt the initial data used by GR03b for easier comparison between their results and ours. A long 
integration yields a Lyapunov exponent i ~ 0.3 yr~^, in agreement with that found by GR and 
other authors (see Figure 9). 

The variations of J and ■0 ^or the satellite system (plotted in Figure 10) are reminiscent of 
those of the analogous variables pertaining to the parametric pendulum as illustrated in Figure 
4. Similar patterns of switching between libration and circulation are exhibited by each system. 
Fractional variations of the action are larger in the satelUte system because of its lower adiabaticity. 



4.1. Where separatrix crossings occur 

GROSab and subsequent authors state that jumps in the mean motions occur near apse an- 
tialignment. Apse antialignment corresponds to the time of maximum amplitude of the potential. 

We sec in Figure 10 that separatrix crossings currently occur in that region. A similar behavior is 
seen in Figure 8a, in which the action is close to its maximum value in the chaotic zone. Jumps 
in mean motion correspond to intervals spent in libration amidst stretches of circulation during 
which il^ increases monotonically. Changes in the slope, ^, from one episode of circulation to the 
next reflect changes in the action at the intervening separatrix crossings. The jumps are not due 
to chaos; rather they are due to the separatrix crossings, which also give rise to chaos. Changes 
in the circulation rate can be predicted, but only with an error that increases at every separatrix 
crossing. 

A chaotic system explores the entire chaotic zone in phase space. We expect the Prometheus- 
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Pandora system to sample states that are unlike its current state in which the separatrix crossings 

occur near apse antialignment as seen in Figure 8a. Averaged over time, separatrix crossings 
should not be restricted to a narrow range of processional phase. In particular, states in which 
the separatrix crossings take place near apse alignment, such as shown in Figure 8b, must also 
be sampled. Results from our integrations displayed in Figure 11 confirm this expectation; the 
system experiences both states in which it is mostly in circulation and those in which it is mostly 
in libration. Separatrix crossings, accompanied by sudden changes in the behavior of tpit), happen 
at all precessional cycles. Switching between regimes of behavior does not require a "kick" as such; 
the requirement is that the potential change so that the system action moves to the other side of 
the separatrix action. This can happen as easily near apse alignment as near apse antialignment. 

A long integration confirms that the action wanders throughout the chaotic zone (Figure 12, 
cf. also Figure 6). Eqs. (17) and (18) predict that it takes ~ 170 years for the action to diffuse all 
the way across the zone. Due to the weak adiabaticity in this case, a single large jump can span a 
large fraction of the total range of action. Flipping between moderately high and moderately low 
action states can therefore occur in a much shorter time, as seen in Figure 12. 

Within measurement errors, a variety of orbits fit the best orbital data for Prometheus and 
Pandora. Renner et al. (2005) find that after two separatrix crossings (their 2004 point), there 
is an uncertainty of 0.2° in the longitude of each satellite.^ Starting from this level of positional 
uncertainty and holding all other orbital elements fixed, we integrate forward in time a set of 3 
trajectories spanning the error range. The results are shown in Figure 13, in which the action for 
each trajectory is plotted against the separatrix action. In only 15 years from the start of our 
integration, one of the trajectories has already transitioned from a state of high to low action (due 
to experiencing a large jump). This illustrates that only a small uncertainty in positions can rapidly 
amplify to give a qualitative change in system behavior. 



4.2. Limitations of the analogy 

The Prometheus-Pandora system is not strongly adiabatic, since only about 3 periods of os- 
cillation take place during the modulation of the potential. Because the amplitude of the potential 
varies substantially over the period of motion, we expect our frozen potential assumption to lead 
to uncertainties of order half a period of the motion in the positioning of the separatrix crossings. 
However, we have seen that the crossing positions are adequately located by this method. 

A more serious problem is the fact that we do not fulfill the condition (f) <^\V{t)\ throughout 
all of the precessional cycle. The dashed line plotted in Figure 1 shows the magnitude of the cf) term 
in comparison with \V{t)\ = IdV/dtpl- For about one third of the time (around apse alignment), 
we do not fulfill this condition, and the pendulum equation does not describe the system well at 



*This is about the same as the size of a single HST error box at Saturn. 
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all. During these times, the system is essentially governed by ~ —(f), so * appears to circulate. 
Perhaps then we cannot define a separatrix crossing near apse alignment. We do find some evidence 
(not detailed in this paper) that greater jumps in the separation of trajectories occur near apse 
antialignment than near apse alignment. Undeniably however, changes in the evolution of ■i/'lO 
occur throughout the precessional cycle (Figure 11), and this statement is independent of the 
pendulum analogy, since it is based on integrations of the full equation of motion, Eq. (5). The 
behavior of the orbits between apse antialignments certainly cannot always be fitted using freely 
precessing ellipses. 

5. Summary and conclusions 

The equations of motion of the Prometheus-Pandora resonant system are well represented by 
those of a parametric pendulum. Despite the low degree of adiabaticity of the system, the same 
regimes of behavior exist as found in a more adiabatic model system. Using this analogy, we have 
explained the nature of the "kinks" in mean longitude seen in both observations and simulations as 
short episodes of libration between long stretches of circulation. We caution that the system will 
not continue indefinitely to display periods of drift with intervening kinks: in as little as 15 years, 
the evolution of the mean motions could be drastically different. Sudden changes in behavior (and 
increased uncertainties in predictions) can happen far from times of apse antialignment. 
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Quantity 


Prometheus 


Pandora 


Symbols 


unprimed 


primed 


m/M 


5.80 X 10"^° 


3.43 X 10-1° 


xn 


188.53815 


82.14727 


n (° s-i) 


6.797331 X 10-3 


6.629506 X 10-^ 


e 


2.29 X 10-3 


4.37 X 10-3 


wn 


212.85385 


68.22910 




3.1911 X 10-5 


3.0082 X 10-5 



Table 1: Properties of Prometheus and Pandora and of their orbits, from GR03b 
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Tables (continued): 



Resonance q Period (yr) = —2'K/tpq Coefficient Cq 

1 1.078 -1.08 X 10-3 

2 1.303 6.26 X 10-^ 

3 1.648 -1.21 X 10-2 

4 2.239 7.82 x 10-^ 



Table 2: Properties of the four components of the 121:118 mean motion resonance, from GR03b 
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Properties of the potential in which the Prometheus-Pandora system moves: (a) 
soHd hne-the amphtude of the potential as a function of time, dashed Ime-cf); (b) 
Drift phase ^ as a function of time; (c) Rate of change of the drift phase. This is the 
drift rate of the hills in the potential. In most of this study we ignore the variation 
of 4>, since it is localized to a small fraction of the period 16 

The potential V{tp) in which the simplified parametric pendulum moves, for a = 0.2. 
Shaded area shows range over which potential varies with time 16 

For the model parametric pendulum system with a = 0.2, k = 0.01: (a) Action as 
a function of time; thick line is separatrix action and thin line is system action, (b) 
System period as a function of time, with vertical lines at locations of separatrix 
crossings, where the period becomes large 17 

For the same system as in Fig. 3: (a) Zoom-in on system action, showing jumps 
in action at separatrix crossings (marked by vertical lines), (b) Evolution of tp{t), 
showing transitions from libration to circulation and vice versa at separatrix crossings. 17 

For the same system as in Fig. 3, we illustrate the separation of neighboring tra- 
jectories in phase space, (a) A plot of the phase space separation as a function of 
time shows the jumps in both the separation and the rate of separation that occur at 
separatrix crossings, (b) A longer integration gives the system Lyapunov exponent 
i as the asymptotic value of ln(separation)/time at late time, which from the plot is 
e ~ 0.002 18 

Using a less adiabatic model parametric pendulum than in Fig. 3 (k = 0.1, a = 0.2), 

we plot the value of the system action over about 5000 separatrix crossings. The 
action diffuses more rapidly than that in Fig. 3, allowing us to more efficiently 
illustrate wandering across the entire chaotic zone 18 

The system action diffuses through the chaotic zone due to quasi-random jumps in 
action at separatrix crossings. The jump size is given by Eq. 16, and assuming a 
uniform distribution of sin(V's/2), we can construct a probability distribution P(AJ) 
for jump size. The contribution to the diffusivity from a given jump size AJ is 
proportional to (AJ)^P(AJ). This quantity is plotted in the figure, normalizing 
jump size to y = iA/A. Jumps of more than 2 times this size contribute about half 
of the diffusivity while only contributing about 4 % by number 19 

Evolution of the pendulum in Fig. 6, showing the behavior of when the action is 
close to the (a) top and (b) bottom of the chaotic zone 19 
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9 Determination of the Lyapunov exponent for the Prometheus-Pandora system: this 
plot is analogous to Fig. 5b for the model parametric pendulum, in which the 
Lyapunov exponent is given by the late-time value of ln(separation)/time 20 

10 Action and ip for the Prometheus-Pandora system for the first 12 years of integration, 
showing changes in the regime of motion at scparatrix crossings, as indicated by 
vertical lines. The dashed line denotes the separatrix action. The unit of action for 

the Prometheus-Pandora system is rad^ yr~^ throughout 20 

11 Examples of high, medium, and low action states for the Prometheus-Pandora sys- 
tem. Vertical lines indicate times of apse antialignment, not separatrix crossings, 
which themselves occur at different precessional phases in these three cases. Dashed 
lines are separatrix action 21 

12 Wandering of the action of the Prometheus-Pandora system: analogous to Fig. 6 for 

a model pendulum, this clearly delineates the chaotic zone 22 

13 We start with three sets of initial conditions of Prometheus and Pandora correspond- 
ing to 0.2° uncertainty in the longitude of each satellite, and plot the action (along 
with the separatrix action, thick line) for each of these trajectories. In as little as 15 
years, one of these trajectories has transitioned from a high to a low action state. . . 22 
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Fig. 1. — Properties of the potential in which the Prometheus-Pandora system moves: (a) sohd 
hne-the amphtude of the potential as a function of time, dashed line-cp; (b) Drift phase 4> as a 
function of time; (c) Rate of change of the drift phase. This is the drift rate of the hills in the 
potential. In most of this study we ignore the variation of (p, since it is localized to a small fraction 
of the period. 
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Fig. 2. — The potential V{ip) in which the simplified parametric pendulum moves, for a = 0.2. 
Shaded area shows range over which potential varies with time. 
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Fig. 3. — For the model parametric pendulum system with a = 0.2, k = 0.01: (a) Action as a 
function of time; thick line is separatrix action and thin line is system action, (b) System period as 
a function of time, with vertical lines at locations of separatrix crossings, where the period becomes 
large. 
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Fig. 4. — For the same system as in Fig. 3: (a) Zoom-in on system action, showing jumps in action 
at separatrix crossings (marked by vertical lines), (b) Evolution of ip{t), showing transitions from 
libration to circulation and vice versa at separatrix crossings. 
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Fig. 5. — For the same system as in Fig. 3, we illustrate the separation of neighboring trajectories 
in phase space, (a) A plot of the phase space separation as a function of time shows the jumps 
in both the separation and the rate of separation that occur at separatrix crossings, (b) A longer 
integration gives the system Lyapunov exponent i as the asymptotic value of ln(separation)/time 
at late time, which from the plot is £ ~ 0.002. 
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Fig. 6. — Using a less adiabatic model parametric pendulum than in Fig. 3 (k = 0.1, a = 0.2), we 
plot the value of the system action over about 5000 separatrix crossings. The action diffuses more 
rapidly than that in Fig. 3, allowing us to more efficiently illustrate wandering across the entire 
chaotic zone. 
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Fig. 7. — The system action diffuses through the chaotic zone due to quasi-random jumps in action 
at separatrix crossings. The jump size is given by Eq. 16, and assuming a uniform distribution of 
sin(^s/2), we can construct a probabihty distribution P(AJ) for jump size. The contribution to 
the diffusivity from a given jump size A J is proportional to (AJ)^P(AJ). This quantity is pfotted 
in the figure, normaUzing jump size to y = A^Aj A. Jumps of more than 2 times this size contribute 
about half of the diffusivity while only contributing about 4 % by number. 
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Fig. 8. — Evolution of the pendulum in Fig. 6, showing the behavior of when the action is close 
to the (a) top and (b) bottom of the chaotic zone. 
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Fig. 9. — Determination of the Lyapunov exponent for the Prometheus-Pandora system: this plot 
is analogous to Fig. 5b for the model parametric pendulum, in which the Lyapunov exponent is 
given by the late-time value of ln(separation)/time. 
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Fig. 10. — Action and tp for the Prometheus-Pandora system for the first 12 years of integration, 
showing changes in the regime of motion at separatrix crossings, as indicated by vertical lines. The 
dashed line denotes the separatrix action. The unit of action for the Prometheus-Pandora system 
is rad^ yr~^ throughout. 
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Fig. 11. — Examples of high, medium, and low action states for the Prometheus-Pandora system. 
Vertical hues indicate times of apse antialignment, not separatrix crossings, which themselves occur 
at different precessional phases in these three cases. Dashed lines are separatrix action. 
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Fig. 13. — We start with three sets of initial conditions of Prometheus and Pandora corresponding 
to 0.2° uncertainty in the longitude of each satellite, and plot the action (along with the separatrix 
action, thick line) for each of these trajectories. In as little as 15 years, one of these trajectories 
has transitioned from a high to a low action state. 



